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ABSTRACT 

Using a population synthesis technique, we have calculated detailed models of the present-day field 
population of objects that have resulted from the merger of a giant primary and a main-sequence or 
brown dwarf secondary during common-envelope evolution. We used a grid of 116 stellar and 32 low- 
mass/brown dwarf models, a crude model of the merger process, and followed the angular momentum 
evolution of the binary orbit and the primary's rotation prior to merger, as well as the merged object's 
rotation after the merger. We find that present-day merged objects that are observable as giant stars 
or core-helium burning stars in our model population constitute between 0.24% and 0.33% of the initial 
population of ZAMS binaries, depending upon the input parameters chosen. The median projected 
rotational velocity of these merged objects is ~ 16 km sec -1 , an order of magnitude higher than the 
median projected rotational velocity in a model population of normal single stars calculated using the 
same stellar models and initial mass function. The masses of the merged objects are typically less than 
~ 2M Q , with a median mass of 1.28 Mq, which is slightly more than, but not significantly different 
from, their normal single star counterparts. The luminosities in our merged object population range 
from <~ 10 — 100 L@, with a strong peak in the luminosity distribution at ~ 60 L©, since the majority 
of the merged objects (57%) lie on the horizontal branch at the present epoch. The results of our 
population synthesis study are discussed in terms of possible observational counterparts either directly 
involving the high rotational velocity of the merger product or indirectly, via the effect of rotation on 
envelope abundances and on the amount and distribution of circumstellar matter. 
Subject headings: binaries: close — circumstellar material — stars: horizontal branch — stars: rotation 



1. INTRODUCTION 

There is a general consensus that the common-envelope 
(CE) phase plays an essential role in the evolution of 
many close binary stars. During CE evolution, one star 
becomes engulfed in the envelope of its companion. The 
orbit of the binary d ecays due to the acti on of gravi- 
tational torques (e.g., iRicker fe Taam 2008J) between it 
and the non-co-rotating envelope. CE evolution involv- 
ing stars in various evolutionary states can occur, but the 
ultimate fate of such evolution, regardless of the states 
of the component stars, is either a merger or a stable bi- 
nary system (see llben fe Liviolll993l iTaam fe SandquisH 
120001 for reviews). 

There has been a renewed interest in CE mergers, 
specifically in massive binaries where merger of the 
cores or of a core and in-spiraling compact object has 
been suggested as a mechanism for producing hyper- 
novae and gamma-ray bursts (e.g.. IFrver fe Wooslevl 
[19981 Middleditch 2004] IFrver fc Hegerll2005l) . In this 
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paper, we investigate the evolutionary consequences of 
CE mergers involving lower-mass binaries consisting of 
a low- to intermediate-mass giant star and a stellar or 
substellar companion. Research on CE evolution in- 
volving low-mass binaries has focused almost exclusively 
on systems that survive CE evolution. In part, this is 
due to the identification of observational counterparts to 
such surviving systems [e.g., cataclysmic variables, dou- 
ble degenerate white dwarfs) and the lack of such obvi- 
ous counterparts to mergers. However, detailed hydrody- 
namical calculations of CE evolution, incomplete as they 
are, i ndicate that mergers should occur in such systems 
(e.ff.,ISandquist et al.ll998llSandauist. Taam. fc Burkertl 
I2000T ). Furthermore, population synthesis of close bi- 
nary stars that involve a single phase of CE evolution 
between a low- to intermediate-mass giant star and a stel- 
lar or substellar companion predict that mergers should 
be as fr equent as post-CE bina ries in the disk of our 
Galaxy ([Politano fc Weileil 120071) . Such studies suggest 
that there may exist a significant number of merged ob- 
jects in our Galaxy that have yet to be recognized as such 
observationally. 

Only a handful of studies of the consequences of CE 
mergers in low-mas s binaries exist in the literature. 
IPodsia dlowski (200lJ) provides an overview of some of the 
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physical consi derations involved in mo deling CE merg- 
ers in genera l. iSiess fc Livid (11999allbT). iSoker fc Harpazi 
((20001 12001 . and ICarlberg. Maiewski. fc Arrasl (|2009h 
investigated the spiral-in and evaporation of planets or 
brown dwarfs (BDs) inside the envelope of a giant star. 
These authors found that the result of such a merger 
is a rapidly rotating single star, and they suggest that 
such mergers could explain certain abundance anoma- 
lies, such as enhanced Li, observed in the envelopes of 
a small fraction of giant stars, as well as the hetero- 
geneity found among horizontal branch (HB) stars. Blue 
stragglers, single subdwarf B (sdB) or subdwarf O (sdO) 
stars, and R Coronac Borcalis stars, to name a few, have 
been suggested as possible observational counterparts to 
low-mass CE mergers (e.g.. iTutukov fc Yunglesof] 120051 : 
ISoker fc Harpad l2007t IPolitano et al.1 120081 ) However, 
there have been only two population synthesis studies 
that explored the consequences o f low- mass CE merg- 
ers: llzzard. Jefferv. fc Lattanziol (|2007l ) , who modeled 
the po pulation of early- type R stars, and IPolitano et all 
(|2008|) . who modeled the population of single sdB stars. 
In this paper, we extend our previous work 
(|Politano et al.ll2008[ ) and present population synthesis 
calculations of the present-day population of objects in 
the disk of our Galaxy that have resulted from merg- 
ers during a single phase of CE evolution involving a 
giant primary and a stellar or substcllar secondary To 
the best of our knowledge, these are the first popula- 
tion synthesis calculations that comprehensively explore 
the consequences of mergers that occur during CE evo- 
lution involving low-mass binaries. In §[2J we describe 
the method used to model this population and discuss 
the assumptions underlying our study. In §02 the re- 
sults of our calculations are presented and compared to 
a corresponding present-day population of normal single 
stars calculated using the same grid of stellar-evolution 
models. In §[U we discuss the observational signatures of 
CE mergers and speculate about potential observational 
counterparts to CE mergers. Finally, we summarize and 
conclude in the last section (§[5]). 

2. METHOD 
2.1. Population synthesis code 

We use the Monte- Carlo population synthesis code de- 
velop ed by Politano (|Politanol 119961: IPolitano fc Weilerl 
|2007f ) with two major modifications : (1) the analyt ic fits 
that had been used previously (see lPolitanolll996T ) have 
been replaced by numerical tables of 116 up-to-date stel- 
lar models ranging in mass from 0.5 to 10.0 Mq in incre- 
ments of 0.1 Mq and from 10.5 to 20 Mq in increments 
of 0.5 Mq and (2) tidal effects which act to synchronize 
the rotation of the primary with the orbital motion of 
the system have been included. 

The updated stellar models were calculated using a ver- 
sion of the bi nary stellar evo l ution c ode e-jj developed 
by E ggleton (jEggletonl I1971L Il972t lYakut fc Eggletonl 
|2005|. and reference s therein) and updated as described 
in iPols et al.l (j!995t ). Convective mixing is modelled by 
a diffusion equation, using a mixing-length to pressure 
scale-height ratio l/H p = 2.0. Convective overshooting 

4 The current version of ev is obtainable on request from 
eggletonlOllnl . gov, along with data files and a user manual. 



in the core is taken into account on the main sequence 
for stars with M > 1.2 Mq and on the horizontal branch 
for stars of all masses. We use an overshooting param- 
eter 6 ov = 0.12, which corresponds to an overshooting 
length of about 0.3 H p . The code cannot evolve a model 
through the helium flash, but this is resolved by auto- 
matically replacing the model at the moment of degen- 
erate helium ignition by a tailored model with the same 
total mass and core mass in which helium was just ig- 
nited non-degenerately. For stars with M > 2.1 Mq, 
helium ignites non-degenerately and this intervention is 
not needed. The helium-core mass M c is defined as 
the mass coordinate at which the hydrogen abundance 
drops below 10%. We compute the binding energy of the 
hydrogen-rich envelope, -E^ind, by integrating the gravi- 
tational and internal energies over the mass coordinate 
of the model, from the core-envelope boundary to the 
surface of the star. In the internal energy we include 
the thermal energy, but not the recombination energy. 
M ore details regardin g thes e assumptions are provided 
in Ivan der Sluvs et al.l (|2006() . The initial composition of 
our model stars is X = 0.70, Y = 0.28 and Z = 0.02. 
Mass loss via stellar winds is incorporated in the mod- 
els using a pre scription that was originally inspired by 
iReimersI ((19751) : 
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where we use r\ = 0-2 in this study, and the prescription 
by Ide Jager et al.l (J1988D . which is negligible except for 
the most massive stars in our grid. The upper prescrip- 
tion in Eq.Q]is used in the majority of our models, except 
where the envelope binding energy is small (in absolute 
value), e.g., for stars near the tip of the asymptotic giant 
branch (AGB). 

As a second modification, we now assume synchronism 
of the primary's rotation with the orbital motion when- 
ever the primary is on either giant branch. In these cases, 
after each time step, the new total angular momentum 
of the system is calculated, accounting for losses due to 
stellar winds. Synchronism is then enforced by redis- 
tributing this total angular momentum, Jtot, such that 
the rotational angular velocity of the primary is equal to 
the orbital angular velocity of the binary. The new or- 
bital separation is calculated as the root of the following 
equation: 



bz 3 + c = 0. 



(2) 



In this equation, z = (a/R) 1 ' 2 , where a is the orbital 
separation and R is the primary's radius, and b and c 
are dimensionless constants given by b = Jtot/[G/J, 2 (M p + 
M S )R} 1 / 2 and c = k 2 (l + M p /M s ), where G is the univer- 
sal gravitation constant, M p and M s are the primary and 
secondary masses, respectively, fi is the reduced mass of 
the binary, and k is the radius of gyration of the primary. 
In all cases, the rotational angular momentum of the sec- 
ondary is neglected and a circular orbit and rigid-body 
rotation for the primary are assumed. 

In our population synthesis calculations, we begin with 
10 7 zero-age main sequence (ZAMS) binaries. We assume 
that these binaries arc distributed over primary mass 
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(M p ) according to a iMiller fc Scald (|1979f ) initial-mass 
function (I MF), over orbital period uniformly in logP 
(|Abtl I1983T ) an d over mass ratios uniformly in q (i.e., 
q(q) da = 1 da. iDuauennoy fe Mavorlll99H iMazeh et al.1 
11992: IGoldberg et al.l 120031) . where a = M s /M p . Given 
the large observational uncertainties in the distribution 
of mass ratios in ZAMS binaries, we also investigate 
the dependence of our results on the assumed choice for 
g(q) (see § 12. 4j) . We adopt a minimum primary mass 
of 0.95 Mq, a maximum primary mass of 10 M Q and 
a minimum secondary mass of 0.013 M@ in our calcu- 
lations. For secondaries with masses less than 0.5 M©, 
including substellar secondari es, we use detailed stellar 
models from the Lyon g r oup (IChabrier fc Baraffei Il997t 
iBaraffe et afl[T99l [20031: IChabrier et al.H2000D . For sec- 
ondaries with masses greater than or equal to 0.5 M@, 
we use the same stellar models as for the primary. We 
assume that the CE phase is so rapid that its dura- 
tion is negligible compared to the other time scales (e.g., 
Ilben fe Livid [llMlTaarn fe Sandauistl 120001) . the age of 
the Galaxy is 10 10 yrs (representative of the thin-disk 
population) , and the star formation rate throughout the 
Galaxy's history has remained constant. Lastly, we do 
not take into account angular momentum loss due to 
magnetic braking. 

2.2. Merger scenario 

To model the population of mergers between giant pri- 
maries and less massive companions, we adopt the follow- 
ing scenario. A given ZAMS binary is evolved until the 
primary reaches the base of the red-giant branch (RGB), 
at which point synchronism between the primary's rota- 
tion and the orbit is assumed. As the primary ascends 
the RGB, two conditions are checked at the end of each 
time step: (1) has the primary filled its Roche lobe or 
(2) has a tidal instability developed (see below). If the 
Roche lobe has been filled, we check to sec if the mass 
transfer will be uns table by using the criterion given in 
IHurlev et all (120021) (their Eq. 57): 
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where M p and M c are the total mass and core mass of 
the primary when it fills its Roche lobe and x = 0.30 if 
M c < 0.50 M Q or x = 0.33 if M c > 0.50 M . If q < 
q cl -it, the mass transfer will be unstable and a CE phase 
will ensue. If q > q C mt, the mass transfer will be stable 
and the system will avoid a CE phase. The evolution of 
such systems is no longer followed and these systems are 
removed from the population^. If a tidal instability has 
developed, we assume a CE phase will inevitably ensue. 
If the primary does not fill its Roche lobe or if a tidal 
instability does not develop on the RGB, the primary 
is evolved through the core He-burning phase and then 
resynchronized at the base of the AGB. As the primary 
ascends the AGB, the same two conditions as before are 
checked at the end of each time step and the same pro- 
cedure is followed. If neither condition is satisfied when 

5 We recognize that such systems may potentially undergo a 
CE phase later in their evolution (e.g., if unstable mass transfer 
is initiated when the secondary fills its Roche lobe). However, our 
code does not treat phases of stable mass transfer. 



the primary reaches the tip of the AGB, the system is 
removed from the population. 

To determine if merger occurs within the CE, sim- 
ple energy considerations are us ed to relate the pre- 
and p ost-CE orbital separations (jTutukov fc Yunglesonl 
119791 ). Although the exact prescription varies somewhat 
depending upon the author, a typical expression used 
in population sy n thesis calculatio n s is gi v en below (e.g., 
Iben fc Tutukovl IT9851 : iPolitand H99a iWillems et all 
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where a; and at are the pre- and post-CE orbital sepa- 
rations, respectively, -EWd.i is the binding energy of the 
primary's envelope at the onset of the CE phase, and 
«ce is a parameter that measures the efficiency with 
which orbital energy is transferred to the CE. The pa- 
rameter «ce embodies a major uncertainty in this sim- 
plified prescription. We have chosen acE = 1 as our 
standard model and we investigate the dependence of 
our results on acE (see § 12. 4j) . The term i?bind,i includes 
both the gravitational binding energy and thermal in- 
ternal energy of primary's envelope, and is computed 
directly from the detailed stellar-structure models. We 
note that this eliminates the need to parametrize the en- 
velope's binding energy using a dimensionlcss constant, 
A (i.e., S b i„d,i = GM C [M C + M c ]/XR h where M c and R { 
are the envelope mass and radius of the primary at the 
onset of the CE phase). Merger is assumed to have oc- 
curred if the radius of the secondary is larger than its 
Roche-lobe radius at the end of the CE phase. 

We also consider mergers that result from a tidal insta- 
bility. If the moment of inertia of the primary exceeds 
one-third of the moment of inertia of the binary orbit 
before the primary fills its Roche lobe on the RGB or 
AGB, Eq.[2] has no real roots. Physically, this implies 
that there will no longer be sufficient angular momen- 
tum within the orbit to k eep the envelope rotating syn- 
chronously (|Darwinlll879r) and the orbit will continue to 
shrink until mass transfer starts, after which either the 
envelope is ejected or a merger occurs. 

2.2.1. Treatment of the merger process 

Available detailed hydrodynamical models of stellar 
mergers have focused on collisions between stars and 
are more appropriate for mergers that occur as a result 
of stellar interactions within dense stellar env ironments, 
such as at the centers of globular c lusters (e.o.. lSills et al.1 
119971 I200H lLombardiet al.ll2002| ). In the absence of de- 
tailed models that arc applicable to stellar mergers dur- 
ing CE evolution, a very simplified model for the merger 
process is adopted. The key assumptions underlying our 
treatment are listed below. 

1. In most cases, the merger of part or all of the 
secondary with the primary leads to an object 
spinning at several times its break-up rotational 
velocity, given by t>b r = (GM/R) 1 / 2 , where M 
and R are the object's mass and radius, respec- 
tively. Consequently, we assume that a rapid mass- 
loss phase occurs when the rotational velocity of 
the object exceeds some critical rotational velocity, 
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Vcrit = /3^br- For our standard model, we assume 
(3 = 1/3, which corresponds to centrifugal forces 
contributing about 10% of the force against grav- 
ity. Our choice is motivated by the fact that we 
have adopted mass loss rates for single stars, in- 
cluding giant stars, that are slowly rotating. At 
higher rotational velocities, it is likely that strong 
magnetic activity would be present, significantly 
enhancing the mass loss rates. The rate of mass loss 
via such a process is not well understood, and we 
have adopted standard rates to avoid introducing 
an additional parameter in this preliminary study. 
To explore the qualitative effect of the choice of /? 
on our results, we calcu late d model populations for 
other values of /3 (see § 12. 4[) . 

2. The merged object is assumed to assimilate only 
as much mass from the binary as is required such 
that it rotates at the critical velocity specified un- 
der 1, and we assume that the remaining mass is 
expelled from the system. We note that the latter 
is a simplifying assumption. In reality, it may be 
that only a portion of this remaining mass actually 
has enough energy to escape from the system, with 
the rest eventually accreted by the merged object 
(likely through an accretion disk). 

3. The material in the merged object is assumed to 
be grossly distributed according to its specific en- 
tropy, with lower-entropy material residing in the 
core and highcr-cntropy material residing in the en- 
velope. Consequently, we assume that the merged 
object has a core mass that is equal to that of the 
original primary, but an envelope mass that is po- 
tentially larger than that of the primary. 

4. Although the merged object is unlikely to be in 
a state of thermal equilibrium initially, we assume 
for definiteness that the radius of the merged ob- 
ject corresponds to that of a star in the model 
grid for the same core mass and total mass after 
it achieves thermal equilibrium. Such an approxi- 
mation may underestimate the amount of angular 
momentum that the remnant may acquire imme- 
diately after the merger process, which could lead 
to an increased level of magnetic activity due to 
dynamo action and an elevated level of mass and 
angular momentum loss (see item 1 above). 

5. In some cases, the merger results in an object 
with a total mass/core mass combination for which 
no corresponding giant-branch model exists in our 
stellar grid. In the majority of these cases, the core 
mass of the merged object is less than the core mass 
at the base of the corresponding giant branch for 
a stellar model with the same total mass. Rather 
than permit an unphysical evolutionary regression, 
such as an AGB primary reverting to a core He- 
burning star as a result of the merger, we choose 
the stellar model in our grid at the base of the cor- 
responding giant branch with the same total mass 
as the merged object. Clearly this is not optimal, 
but the difference in core mass is generally small 
and, at present, we can do no better. 



6. If the CE takes place when the primary is on the 
RGB, the merger may result in a total- mass/core- 
mass combination that corresponds to a core-He 
burning star in our model grid (i.e., since the crit- 
ical core mass for degenerate He ignition decreases 
with increasing total mass, the additional envelope 
material acquired during merger may reduce the 
critical core mass needed for He ignition). Since it 
is physically reasonable for an RGB star to ignite 
helium in the core, even if induced by a merger 
event, we adopt a core He-burning model for the 
merged object. 

2.3. Evolution of the merged object to the present 

epoch 

The subsequent evolution of the merged object, in- 
cluding its rotational velocity, w rot , is followed until the 
present epoch using the same grid of stellar evolution 
models that we used for the primary prior to merger. In 
addition to angular momentum loss due to stellar winds, 
we assume that angular momentum is removed from the 
merged object's envelope via rotationally enhanced mass 
loss whenever w rot exceeds w cr it- Typically, this occurs 
when there is a significant decrease in the merged object's 
radius, such as during contraction to the HB following He 
ignition. If u rot > f cr it during a given time step, we cal- 
culate the amount of angular momentum that needs to 
be lost so that the merged object will rotate exactly at 
t> cr it and we artificially remove the corresponding amount 
of mass from the envelope. Once O.lMg of material is 
lost, it becomes necessary to shift to the next lowest- 
mass model in our grid with the same core mass in order 
to continue to follow the object's evolution. If there is 
no corresponding lower-mass model with the same core 
mass, we follow the procedure described in assumption 
5 in the previous section. However, in some cases this 
procedure proved unsatisfactory and led to unphysical 
behavior. For example, some objects would continue to 
lose mass, moving along the base of the HB until the 
object moved off the grid. To address this, we softened 
our condition for shifting to the next lowest-mass model 
by letting the object evolve for five additional time steps 
before shifting. We found that if the object's rotational 
velocity dropped below critical during these five time 
steps, its rotational velocity stabilized and the evolution 
could generally be followed for the remainder of its evo- 
lution without problem. If the merged object completes 
its evolution as an AGB star before the present epoch is 
reached, we assume it will be a white dwarf at the present 
epoch, and we remove the object from the population. 

2.4. Input parameters 

For our standard model, we have assumed values for 
the CE efficiency parameter and ZAMS mass ratio dis- 
tribution that are common in the literature, namely, 
q;ce = 1 and g(q) dq = 1 dq. In addition, we have cho- 
sen a value of /3 = 1/3 in our standard model, where 
the parameter (3 represents the fraction of the break-up 
velocity for which we assume rotationally-induce mass 
transfer begins (see item 1 of ? !2.2.ip . Given the large un- 
certainties in these parameters, we have calculated model 
populations for other choices of these parameters. We list 
the parameter choices made for each model in Tabic [T] 
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Table 1 

Input parameter choices for the 
models in this study 
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In this section, we present the results of our population 
synthesis calculations. In § 13. 1[ we discuss the general 
characteristics of the present-day population of merged 
objects, with an emphasis on observable quantities such 
as rotational velocity, luminosity, effective temperature, 
etc. All results in § 13.11 are for our standard model (see 
Table [1]). In § 13.21 we discuss the dependence of our 
results on other choices for these three parameters, as 
listed in Table [1] Finally, in § 13.31 we compare our model 
population of merged objects with a standard single-star 
population calculated from the same set of stellar evolu- 
tionary models. 

3.1. Properties of present-day population of merged 

objects 

3.1.1. Relative numbers and sub-populations 

In our population synthesis calculations, we consider a 
population of 10 7 primordial binaries with initial proper- 
ties as described in § 12.11 In our standard model, ~ 71% 
of these binaries never initiate mass transfer, 13% un- 
dergo stable mass transfer and 16% undergo unstable 
mass transfer, which leads to a CE phase. Of those bina- 
ries undergoing CE evolution, 48% survive the CE phase, 
while 52% merge (8.5% of the original population). 82% 
of the CEs are formed when the primary fills its Roche 
lobe, whereas 18% undergo a Darwin instability. Fur- 
thermore, 87% of the mergers occur when the primary 
is on the RGB, while 13% occur when the primary is on 
the AGB. 

In the evolution following merger, the vast majority 
(97%) of the merged objects will have evolved beyond 
the AGB by the present epoch. These merger remnants 
would be observed as white dwarfs today and are not in- 
cluded in our present-day population of merged objects. 
Only 3% of the merged objects (~ 25, 000 stars or 0.25% 
of our original population) would be observed as non- 
degenerate stars today. The properties of these 25,000 
merged objects are listed in Tabled] 

We have divided the total population of 25,000 present- 
day merged objects into three sub-populations: RGB 
stars, HB stars and AGB stars. The properties of each 
of these sub-populations are also listed in Table [2] We 
find that HB stars are the most numerous in our present- 
day population (57%), followed by RGB stars (37%) and 
then AGB stars (6%). This can be understood intuitively 
from the fact that of the merger remnants observed to- 
day, 99.9% were formed in a merger on the RGB, com- 
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Figure 1. Distribution of the rotational velocities in our stan- 
dard model population of present-day merged objects for three sub- 
populations: RGB stars (dashed curve), HB stars (solid curve) and 
AGB stars (dotted curve). The horizontal axis shows the rotational 
velocity as a fraction of the critical rotational velocity, i> cr i t , which 
is assumed to be ^ of the break-up velocity (see text). 

pared to 0.1% for which the merger took place on the 
AGB. In particular, we find that mergers during CE evo- 
lution typically occur in binaries with primary masses of 
^1 — 2 Mq that fill their Roche lobes on the lower half 
of the RGB at orbital periods between roughly 1 and 30 
days. The distribution of mass ratios at the onset of the 
CE phase is almost uniform between and 1. While RGB 
mergers produce remnants that must evolve through the 
rest of the RGB phase, the relatively long-lasting HB 
phase and the AGB phase, AGB mergers produce rem- 
nants that evolve off the AGB quickly and disappear from 
our present-day population. Of all the RGB mergers, 
only 2.4% of the remnants ignite helium in the core and 
become HB stars as an immediate consequence of the 
merger. Most of the remaining RGB mergers will have 
evolved off the RGB by the time of the present epoch, 
and the distribution of stars over the three branches is 
to a large extent determined by the relative lifetimes of 
the stars in each phase. 

3.1.2. Rotational velocity 

In Fig.Q] we show the predicted distribution of present- 
day rotational velocities in our model population of 
merged objects for the three sub-populations: RGB stars 
(red/dashed), HB stars (green/solid) and AGB stars 
(blue/dotted). The rotational velocity is given as a frac- 
tion of an assumed critical velocity for rotational mass 
loss, u cr it (see § 12.2. ip . For our standard model, i; C rit = 
14.6kms- 1 (M/M ) 1 / 2 (i?/lOOi? )- 1 / 2 . 

We find that the distributions for the AGB, RGB and 
HB, respectively, are increasingly peaked towards high 
V ro t- This occurs because most of the stars are on the 
RGB directly after merger and are rotating rapidly. Dur- 
ing the evolution on the remainder of the RGB, the rota- 
tion of the stars will slow down due to their nuclear evo- 
lutionary expansion. When helium ignites in the core, 
the merged objects evolve to the HB, which causes them 
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Table 2 

Dependence of selected model results upon input parameters 



Model 


JV a 


N b 
iVtot 


M c 
(M Q ) 


v sin i d 
(km/s) 


Fraction with 


M re / 
(M ) 


M roj g 
M bin 


AA/ mrg h 




Drot <0.1v crit 


V rot = f crit 






RGB stars 


standard 


9301 


0.37 


1.20 


18.4 


(0.001) 


0.0319 


0.63 


0.34 


0.00 


oce = 0.5 


9534 


0.34 


1.20 


17.7 


(0.0009) 


0.0343 


0.64 


0.34 


0.00 


«CE =0.1 


9634 


0.29 


1.20 


17.4 


(0.0009) 


0.0363 


0.64 


0.34 


0.00 


aiq) = i 


9044 


0.36 


1.20 


18.2 


(0.0009) 


0.0368 


0.89 


0.42 


0.00 


g(q) = q~ ' 9 


8819 


0.35 


1.20 


18.4 


(0.0007) 


0.0251 


0.16 


0.11 


0.00 


fcilt = V hr 


8025 


0.33 


1.30 


48.7 


(0.0006) 


0.0179 


0.45 


0.25 


0.00 


Uorit = 0.1 «br 


9783 


0.38 


1.19 


5.4 


(0.001) 


0.0250 


0.71 


0.37 


0.00 


single stars 


178651 


0.61 


1.20 


1.9 


0.9627 


0.000 








HB stars 


standard 


14305 


0.57 


1.35 


16.1 


(0.0000) 


0.3149 


0.93 


0.40 


0.12 


"CE = 0.5 


16987 


0.60 


1.34 


16.7 


(0.0000) 


0.3469 


0.91 


0.41 


0.12 


«CE = 0.1 


21043 


0.64 


1.23 


17.3 


(0.0000) 


0.4091 


0.91 


0.42 


0.13 


aiq) = 1 


14267 


0.57 


1.35 


16.1 


(0.0000) 


0.3275 


1.25 


0.49 


0.12 


g(q) = q- 09 


14929 


0.59 


1.35 


16.5 


(0.0000) 


0.3264 


0.32 


0.20 


0.11 


Ucrit = «br 


14943 


0.61 


1.58 


50.5 


(0.0004) 


0.2098 


0.74 


0.32 


0.08 


^crit = 0.1 «br 


14217 


0.56 


1.24 


4.5 


(0.0008) 


0.3366 


1.03 


0.44 


0.12 


single stars 


104979 


0.36 


1.58 


3.2 


0.0886 


0.0021 








AGB stars 


standard 


1435 


0.06 


1.34 


6.0 


0.0683 


(0.0007) 


0.94 


0.42 


0.13 


a C E = 0.5 


1748 


0.06 


1.31 


6.3 


0.0732 


(0.0011) 


0.94 


0.42 


0.13 


oce = 0.1 


2205 


0.07 


1.23 


6.4 


0.0757 


(0.0018) 


0.92 


0.43 


0.14 


aii) = q 


1542 


0.06 


1.34 


6.4 


0.0746 


(0.0013) 


1.27 


0.49 


0.13 


g(q) = g-0-8 


1595 


0.06 


1.34 


6.0 


0.0678 


(0.0031) 


0.36 


0.23 


0.12 


"crlt = u br 


1446 


0.06 


1.56 


18.7 


0.0788 


(0.0007) 


0.74 


0.33 


0.09 


^crit = 0.1 «br 


1490 


0.06 


1.23 


1.8 


0.0859 


(0.0007) 


1.09 


0.46 


0.14 


single stars 


10487 


0.04 


1.45 


1.3 


0.5657 


(0.0000) 








Total population 


standard 


25041 


1.00 


1.28 


16.2 


0.0043 


0.1918 


0.81 


0.38 


0.07 


a C E = 0.5 


28269 


1.00 


1.23 


16.2 


0.0048 


0.2201 


0.81 


0.38 


0.08 


«CE = 0.1 


32882 


1.00 


1.23 


16.5 


0.0054 


0.2726 


0.83 


0.40 


0.10 


aiq) = q 


24853 


1.00 


1.29 


16.0 


0.0049 


0.2015 


1.10 


0.46 


0.07 


aiq) = q~ ' 9 


25343 


1.00 


1.29 


16.3 


0.0045 


0.2012 


0.28 


0.18 


0.07 


*W = «br 


24414 


1.00 


1.46 


47.7 


0.0051 


0.1343 


0.63 


0.30 


0.02 


^crit = 0.1 V bl 


25490 


1.00 


1.20 


1.6 


0.0058 


0.1974 


0.90 


0.41 


0.08 


single stars 


294117 


1.00 


1.23 


2.3 


0.6366 


0.0008 









Note. — Entries enclosed by parentheses indicate the number of objects is so small as to not be statistically 
significant. 

a total number of merged objects in present-day population 
b fraction of present-day population in this evolutionary state 
c median mass of merged objects in present-day population 

d median projected rotational velocity of merged objects in present-day population 

c fraction of merged objects in present -day population rotating at one-tenth of the critical velocity or less and at 
the critical velocity (defined in § 12,2.11 1 

median amount of mass that is not accreted by the remnant during merger 
s same as f, except given as a fraction of the total mass of the binary at the onset of the CE phase 
h median of the relative difference between the mass immediately after the merger event, M mrg j, and the mass 

at the present epoch, Af mrgi pE: — "" g .' f mrg ' PE 



to shrink and increase their rotational velocities. Many 
of these stars will reach critical rotation while shrinking, 
thus forcing them to lose mass and angular momentum as 
described in § 12.31 These stars will arrive on the HB spin- 
ning critically, causing the peak at log^rot/fcrit) ~ 0. 
Since HB stars do not expand dramatically during the 
core helium burning phase, many will be found at or 
close to this peak. When these stars continue their evo- 
lution on the AGB, they will be spinning less rapidly 
than they were on the RGB at a similar radius, because 



of angular momentum loss during the intervening evolu- 
tion. In fact, of the three sub-populations, there exists 
a relatively sharp lower cutoff in rotational velocities at 
0.1w cr it for both the HB and RGB stars, whereas the 
AGB stars show a smooth decline. For merged objects 
that are AGB stars at the present epoch, we find that 
98% were formed through a merger on the RGB, while 
only 2% actually merged on the AGB. 

To facilitate comparison with observations, Figure [5^, 
shows the same distributions as Fig.Q] but with the 
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velocity converted to a physical equatorial velocity in 
kms -1 . In addition, in Fig. [2}) we have computed the 
projected rotational velocity, vsini, as could be inferred 
from spectroscopic observations. We used a uniform dis- 
tribution in sin i to generate a random inclination for 
each of our merger remnants. Whereas most late- type 
giant stars are observed t o rotate with projected veloc- 
ities of a few kms -1 (e.q-. lCariberg. M aicws ki. fc Arra"sl 
|2009( ) — indeed our comparison model of a single-star 
population (see § 13. 3p has a median wsinz of 2.3 kms -1 
— we find that the median value for v sin i of the total 
merged population in Fig. [2}) is 16.2 kms -1 . 

In Figure [3l we show an HR diagram with the distribu- 
tion of the medians of the projected rotational velocities 
in each bin. The highest velocity stars are found at the 
bottom of the RGB, with physical and projected veloc- 
ities in excess of 100 kms -1 and 80 kms -1 , respectively. 
Stars on the RGB around 5000 K and on the HB still have 
projected velocities between 10 and 60 kms -1 , whereas 
the merger products in the upper regions of the HR di- 
agram have rotational velocities that are more common 
for giants (a few kms -1 ), although still faster than for 
single stars in that region of the HR diagram. 

3.1.3. Total mass 

In Figure 31 we show the predicted distribution of 
present-day total masses in our model population of 
merged objects for the three sub-populations: RGB stars 
(red/dashed), HB stars (green/solid) and AGB stars 
(blue/dotted) . All three distributions peak at low masses 
and monotonically decrease towards higher mass, as may 
be expected. However, while the AGB distribution looks 
much like a scaled-down version of the HB distribution, 
the distribution for merged objects presently on the RGB 
is much narrower, with a higher low-mass cut-off and 
a steeper drop-off towards higher masses. This steep 
drop-off occurs because a low-mass star (M < 2.1 M©) 
with a degenerate helium core spends more time on the 
RGB than a higher- mass star (M > 2.1 M©) with a non- 
degenerate helium core. As a consequence, of all merg- 
ers that occurred while the primary was on the RGB, 
42% of them remain on the RGB at the present epoch if 
M < 2.1 M©, whereas only 1.2% of them are presently 
RGB stars if M > 2.1 M©. As a result of this, merg- 
ers in the latter category end up as rapidly rotating HB 
stars about 2.5 times more frequently and have an 80% 
higher chance to be observed as AGB stars today than 
mergers in the former group. Table [5] shows that, for 
our standard model, the median mass of the total pop- 
ulation is 1.28 M©, with a clearly lower number for the 
RGB sub-population due to this steep drop-off. 

Figure [5] displays a two-dimensional histogram of the 
Mt t-Wrot plane, for the total present-day population of 
merged objects (RGB, HB and AGB). This figure in- 
dicates that most of the merger remnants are low-mass 
(< 2 M©), rapidly rotating stars, with a peak in the dis- 
tribution corresponding to a mass and velocity range of 
M ps 1.0 — 1.3 M Q and v rot /v cr it w 0.75 — 1, respectively. 

3.1.4. Luminosity 

In Fig. [HI the predicted distribution of the present-day 
luminosities in our model population of merged objects is 
illustrated. As in the previous figures, the distribution is 



divided into 3 sub-populations: RGB stars (red/dashed), 
HB stars (green/solid) and AGB stars (blue/dotted). 
The HB star distribution ranges in logL/L© from 1.6 
to ~ 3 and is sharply peaked near logL/L© = 1.8. 
The RGB distribution is broader and flatter, ranging in 
logL/L© from ~ 0.4 to 3.3 and having a broad maxi- 
mum near logL/L© = 1.4. The AGB distribution ranges 
in log L/Lq from ~ 1.9 to 4.1 and has multiple maxima, 
with a strong peak at logL/L© ~ 2.0 and a weaker one 
near log L/Lq ~ 2.4. We note that the strong peak in the 
AGB distribution lies near that of the HB distribution, 
thus contributing to the peak in the overall luminosity 
distribution. 

3.1.5. HR diagram 

In Figure [3 we show binned HR diagrams (HRDs) for 
the total model population of merged objects. Fig. [7k. is a 
two-dimensional histogram, in which the grey scale indi- 
cates the number of merged objects per two-dimensional 
bin (pixel). If we draw an imaginary line from the upper- 
left to the lower-right corners of the diagram, this line 
intersects the distribution on the giant branches almost 
perpendicularly. When we start in the upper-left corner 
and follow this line through the distribution, we find that 
there is a gradient in the number density of merged ob- 
jects, and that this density increases as we continue to 
follow the line towards the lower-right corner. This gra- 
dient is the result of the IMF; there are more lower-mass 
than higher-mass stars, and the evolutionary tracks of 
lower-mass stars on the giant branches lie to the lower 
right of those of higher-mass stars. In addition, we find 
that the majority of our merged objects sit in a dense 
clump at the HB (see S l3~l~2|) . with T cff « 4500 - 5000 K 
and L rs 60 L©. 

Figure [JJd displays the medians of the velocity distri- 
butions of merged objects in each bin of the HRD. Since 
the standard deviations of these distributions are almost 
always significantly smaller than the medians, this is a 
relevant estimate of the rotational velocity. Whereas the 
number of stars in each bin is a function of mass, the ve- 
locity is clearly dependent on evolutionary phase. Espe- 
cially for the RGB, it is clear that the rotational velocity 
diminishes towards the tip of the giant branch. The HB 
contains fast rotators, as we saw earlier, and even the 
AGB shows a few pixels with a high rotational velocity, 
although Fig. [7b. indicates that these pixels each contain 
about one star. 

3.1.6. Oblateness 

The high rotational velocities of the present-day 
merged objects (see § I3.1.2|) may well give rise to defor- 
mation of the stars. In order to quantify the oblateness of 
the stars, we assumed they can be treated as MacLaurin 
spheroids, for which the relation between angular fre- 
quency to, density p and oblateness e is given by 



y/^Gp 



1- 



,1/2 



(3-2e 2 )sin -1 (e)-— (1- 



1/2 



'(5) 

(jMacLaurinl ll 742J ) . where we use p = p», the mean den- 
sity of the star, and the eccentricity or oblateness is de- 
fined in terms of the polar and equatorial radius (r p and 
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Figure 2. Histograms of the equatorial rotational velocities of the sub-populations of merged objects on the RGB (dashed curve), HB 
(solid curve) and AGB (dotted curve) in our standard model. The left-hand panel (a) shows the distributions of physical velocities in 
km/s, the right-hand panel (b) shows the projected velocities vs'mi as would be inferred from the Doppler shift. 







-= 5000 4000 3000 

Surface temperature of merged objects at present epoch (K) 

Figure 3. Two-dimensional histogram of the projected equatorial 
rotational velocities, vs'mi (in km/s), of the merged objects in our 
standard model, overlaid on an HRD. 



r e , respectively) as 




(6) 



We note that for a constant density, this function re- 
sembles a straight line for the regime of our interest 
(e < 0.85) and we use a simple least-squares method 
to find that the best slope for a line through the origin 
(e = 0,u! = 0) to describe this part of the function (to 
an accuracy of 2.5% or better) is 0.53. This allows us to 
derive an approximation for the oblateness as a function 
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Total mass of merged objects at present epoch (M ) 

Figure 4. Distribution of the total masses in our standard 
model population of present-day merged objects for three sub- 
populations: RGB stars (dashed curve), HB stars (solid curve) 
and AGB stars (dotted curve). 

of the angular frequency and the density 
lu 0.75 w 



0.53 % /2TGp 



e < 0.85 



(7) 



and to compute the oblateness for each of the merger 
products in our standard population. 

Figure [5^, shows the distribution of the oblateness of 
the stars in our present-day populations of RGB, HB 
and AGB stars. We find that for most stars the eccen- 
tricity is larger than 0.1 and should be observable if the 
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Total mass of merged objects at present epoch (M Q ) 

Figure 5. Two-dimensional distribution of the rotational veloci- 
ties and total masses in our standard model population of present- 
day merged objects. The grey scale represents the total number 
of merged objects per two-dimensional bin. As in Fig.fl] the ro- 
tational velocities are given as a fraction of the assumed critical 
rotational velocity. 
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Luminosity of merged objects at present epoch (L Q ) 

Figure 6. Distribution of the luminosities in our standard 
model population of present-day merged objects for three sub- 
populations: RGB stars (dashed curve), HB stars (solid curve) 
and AGB stars (dotted curve). 



star is not too distant (see e.g. IZhao et "all I2009J ). All 
three distributions peak at the maximum value of the 
oblateness (e « 0.5). While the histograms for the RGB 
and HB drop off rapidly towards smaller oblateness, the 
histogram for AGB stars does so less steeply and there 
are a small number of AGB stars with a small eccen- 
tricity (e < 0.01). Figure [S}} displays the distribution 
of eccentricities on an HRD and shows that the most 
deformed stars are at the base of the red-giant branch, 
the HB and the lower AGB, while the most expanded 
giants have shapes that are closer to spherical. The dis- 
tributions of eccentricities follow those of log(w ro t/w C rit) 
in Figs.[T]and[7p. We stress that our model of a MacLau- 



rin spheroid for the shape of these stars has limited va- 
lidity and the fit in Eq.[7] depends on our choice of /3. 
For example, if (3 = 1, this fit, and the assumption of 
MacLaurin spheroids, would yield physically meaningless 
results. Consequently, our use of MacLaurin spheroids 
is meant to provide rough order-of-magnitudc estimates 
and trends, not precise values. 

3.2. Dependence upon input parameters 

In addition to our standard model, we investigated the 
dependence of our results on three uncertain input pa- 
rameters: the CE efficiency parameter, cuce, the initial 
mass ratio distribution in ZAMS binaries, g(q), and the 
assumed critical rotational velocity at which mass is lost, 
Wcrit (see item 1 of § 12.2. 1[) . Uncertainties in the first two 
parameters plague the majority of population synthesis 
calculations involving close binary systems, while uncer- 
tainty in the third parameter is specific to the present 
calculations. We calculated model populations for two 
other choices of each parameter and have listed the pa- 
rameter choices for each model in Table [TJ In all cases, 
only one parameter was varied at a time, keeping the 
other two at their default values. 

Selected results from these non-standard models for the 
three sub-populations discussed in § 13. H and for the total 
population are shown in Table [2] We find that, with a 
few notable exceptions, varying these input parameters 
has relatively little effect on the overall characteristics of 
the present-day population of merged objects. For exam- 
ple, the fraction of the population that is rotating slowly 
(frot < 0.1 i> C rit) is uniformly small for all choices of input 
parameters: in the RGB and HB sub-populations, this 
fraction is either zero or so small as not to be statistically 
significant; in the AGB sub-population, this fraction is 
7 - 9%; and for the total population, this fraction is 
barely significant (< 1%). Further, with the exception 
of the i> cr it = fbr model, the present-day median mass 
in the total population varies from 1.20 M & to 1.29 M Q 
and the variation within a given sub-population is only 
slightly greater. The median mass in the i> cr it = Vb r 
model is systematically greater than the others because 
the merged object can assimilate more angular momen- 
tum (and hence mass) during the merger before it reaches 
its critical rotation velocity (see assumption 2 in § !2.2.ip . 
As expected, the median projected rotational velocity, 
usini, is strongly correlated with our assumed choice 
for w cr it; a larger value for iJ crit allows a higher median 
projected rotational velocity. There is somewhat greater 
variation among the models in the fraction of present-day 
merged objects rotating at v cr - lt . In the total population, 
this fraction varies by a factor of 2, from 0.13 for the 
w C rit = ^br m °dcl to 0.27 for the acE = 0.1 model, com- 
pared with 0.19 for our standard model. The same trend 
is found in the RGB and HB sub-populations. Virtually 
none of the AGB stars are rotating critically, because of 
the large rotational inertia of the envelopes in these stars. 

The largest variation within the present-day popula- 
tion caused by changing the input parameters is found 
in the total number of merged objects. For the total 
population and within a given sub-population, the total 
number of present-day merged objects varies by roughly 
20% to 50%, with the smallest number of objects found in 
the t'crit = «br model population (24,414 stars, or 0.24% 
of the input population) and the largest (32,882 stars, or 



10 



Politano, van der Sluys, Taam & Willcms 





5000 4000 3000 

Surface temperature of merged objects ot present epoch (K) 



5000 4000 3000 

Surface temperature of merged objects ot present epoch (K) 



Figure 7. Two-dimensional distribution of the luminosities and effective temperatures in our standard model population of present-day 
merged objects. In the left-hand panel (a), the grey scale represents the total number of merged objects per two-dimensional bin (pixel). 
The lightest grey scale, which represents a bin with one star, is easily distinguished from the white background. In the right-hand panel 
(b), the grey scale represents the median of the distribution of rotational velocities (as a fraction of w cr it) per pixel. Temperature increases 
to the left on the horizontal axis to facilitate comparison with observational HR diagrams. 
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Figure 8. Distributions of the oblateness e (see Eq.[6jl in our model population of present-day merged objects. Left-hand panel (a): 
histograms for the sub-populations of RGB stars (dashed line), HB stars (solid line) and AGB stars (dotted line). Right-hand panel (b): 
two-dimensional distribution of the oblateness on an HRD. The stars at the base of the RGB are most deformed, whereas very large giants 
are nearly spherical. 



0.33% of the input population) in the «ce = 0.1 model 
(except for RGB stars, where the largest number is found 
in the v cr it =0.1 t>b r model). Reducing the value of q;ce 
results in a less efficient transfer of orbital energy to the 
CE during spiral-in, thereby increasing the number of 
systems that merge during the CE phase. This increase 
in the number of CE mergers results in a larger num- 
ber of present-day merged objects. As one might expect, 
varying t> cr it has very little effect on the total number of 



present-day merged objects. 

Column 8 of Table [5] lists the mass that is not assim- 
ilated from the binary during merger, Af re j, and column 
9 lists this mass as a fraction of the total mass of the 
binary at the beginning of the CE phase. We find that 
varying the choice of cxqe has essentially no effect on the 
amount of mass rejected during merger. Increasing u cr jt 
decreases Af re j , which is as expected since for larger val- 
ues of Ucrit) the merged object can assimilate a greater 
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fraction of the mass of the binary before it starts to spin 
critically. Perhaps surprisingly at first, varying g(q) has 
the strongest effect on the value of M re j. However, this is 
simply because g{q) determines the ZAMS distribution 
of secondary masses. For g(q) = q, equal mass ZAMS bi- 
naries are strongly favored, resulting in systematically 
larger secondary masses than in our standard model. 
Conversely, for g{q) = q~ 09 , extreme mass ratio ZAMS 
binaries are strongly favored, resulting in systematically 
smaller secondary masses than in our standard model. 
Consequently, for g(q) = q, M IC - 3 is larger than in our 
standard model because the secondaries are more mas- 
sive. On the other hand, for g(q) = q~ 09 , M re j is smaller 
than in our standard model because the secondaries are 
less massive. 

The last column in Table [2] lists the fraction of mass 
lost by the merged objects between the moment imme- 
diately after merger and the present epoch. This frac- 
tion is vanishingly small for present-day RGB stars, but 
is ~12% for objects that are currently on the HB and 
AGB, indicating that this post- merger mass loss occurs 
predominantly near the moment of core helium ignition. 
Two factors contribute to this mass loss: (1) normal 
stellar winds near the tip of the RGB prior to helium 
ignition and (2) rotationally-cnhanced mass loss during 
contraction to the HB after helium ignition. A typi- 
cal RGB star spends most of its time on the lower part 
of the RGB without much mass loss. Then, in a rela- 
tively short time, the star evolves to the tip of the RGB, 
where mass loss becomes significant. For example, in the 
1.3 Mq stellar model in our grid, ~ 0.2 M Q is lost on the 
RGB, with 90% of this mass lost in the last 4% of the 
time spent on the RGB. While wind mass loss on the 
RGB is present whether or not merger occurs, we find 
that for a substantial number of objects that merge on 
the RGB, additional mass loss occurs during the star's 
contraction after helium ignition. Especially in stars that 
undergo a helium flash, where the decrease in stellar ra- 
dius is large, the rotational velocity of the merged object 
becomes critical as the star contracts to the HB. A pe- 
riod of enhanced mass loss ensues, until the star has lost 
a sufficient amount of angular momentum for the rota- 
tional velocity to drop below critical. While we do not 
distinguish between these two contributions in column 
10 of Table [2] we note that the post-merger mass loss 
is lower in the v cr i t = i>br model since the larger value 
for w cr it allows the stars to reach larger rotational veloci- 
ties before they start to experience rotationally-enhanced 
mass loss during contraction to the HB. 

3.3. Comparison with single stars 

In this section, we compare our present-day model pop- 
ulation of merged objects to a model population of nor- 
mal present-day single stars. We generated 10 7 single 
ZAMS stars, drawing their initial masses from the same 
IMF that was used for the merged objects, and deter- 
mined the structure and evolution of these stars using 
the same grid of stellar models that was used to gener- 
ate the merged object population (see S I2.1J) . We con- 
sider only those present-day single stars that are on the 
RGB, HB or AGB for this comparison. In Fig.[SJ we 
show the present-day distribution of rotational veloci- 
ties for our single stars in each of these evolutionary 
phases. As before, the rotational velocity for a given 




(v rot /v crJt ) of single stars at present epoch 

Figure 9. Distribution of the rotational velocities in a compari- 
son model of normal present-day single stars constructed from the 
same grid of stellar evolution models used for the merged popu- 
lation. As in the corresponding plot for the merged population 
(Fig.[TJ, the rotational velocities are given as a fraction of the as- 
sumed critical rotational velocity and the distribution is separated 
into three sub-populations: RGB stars (dashed curve), HB stars 
(solid curve) and AGB stars (dotted curve). 



star is shown as a fraction of the critical rotational ve- 
locity for that star, defined as in our standard model 
(u crit = i t>br)- The present-day rotational velocities for 
single stars were calculated using the following expres- 
sion for rotational velocities in ZAMS stars as a function 
of mass (jHurlev et alj|2000l Eq. 107), 



^rot,ZAMS — 



330 M 33 

15.0 + Af 3 - 45 



km sec 



(8) 



conserving angular momentum up to the base of the 
RGB, and then accounting for angular momentum loss 
from the star via stellar winds throughout the remain- 
der of the evolution. Once again, angular momentum 
loss due to magnetic braking was not incorporated in 
the evolution. 

Comparison of Fig. [9] with the corresponding plot for 
the merged objects (Fig.[TJ reveals a striking contrast 
in the rotational velocities for the two populations. 
The peak rotational velocity for each sub-population of 
merged objects is approximately an order of magnitude 
higher than the corresponding peak rotational velocity 
for normal single stars. This difference can be seen even 
more clearly in Fig.QIJl which shows the distribution of 
rotational velocities for the total populations of merged 
objects and single stars. Not only are the peak rotational 
velocities markedly different, but the overlap between 
the two distributions is relatively small. Consequently, 
observational determinations of rotational velocities in 
RGB, HB and AGB stars may be expected to provide an 
important diagnostic in determining which stars are the 
result of a CE merger. 

In Fig. 111! we show the present-day distribution of 
masses for RGB, HB and AGB stars in our single star 
population. Comparison of this plot with the correspond- 
ing plot for the merged objects (Fig. [4} reveals that the 
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(v t /v crit ) of merged objects at present epoch 

Figure 10. Distribution of the rotational velocities in our model 
population of present-day merged objects (solid curve) and in a cor- 
responding present-day population of normal single stars (dashed 
curve) calculated using the same grid of stellar evolution mod- 
els. The solid curve represents the sum of the combined sub- 
populations of RGB, HB and AGB stars shown in Fig. [l] and the 
dashed curve represents the sum of the combined sub-populations 
of RGB, HB and AGB stars shown in Fig.[9] As in those two fig- 
ures, the rotational velocities here arc given as a fraction of the 
assumed critical rotational velocity. 



merged objects have a slightly steeper drop-off towards 
higher masses than the single stars, which is due to 
rotationally-enhanced mass loss from many of the merger 
products during their post-merger evolution. However, 
as can be seen by comparing the median present-day 
masses listed in Table [2] for the total population in our 
standard model and the single-star population, this dif- 
ference is small (0.05 Mq). On the other hand, the rel- 
ative percentages of RGB, HB and AGB stars are quite 
different for the two populations. While the merged ob- 
ject population consists of 37% RGB stars, 57% HB stars 
and 5.7% AGB stars for our standard model, the respec- 
tive percentages for the normal single star population 
are 61%, 36% and 3.6%. The ratios of these percentages 
are 0.61, 1.6 and 1.6, respectively, which suggests that 
this is caused by a shift in numbers between the pre- and 
post-helium-ignition phases. In other words, in a present- 
day population of merged objects on the RGB, HB and 
AGB, the majority of stars have ignited helium, while in 
a corresponding normal single-star population consisting 
of the same three sub-populations, the majority of stars 
have not yet ignited helium. This is understandable since 
for the merged objects, we only consider a star after it 
has merged, which is in 99.9% of the cases somewhere on 
the RGB (see § I3.1.1[) . Hence, its subsequent evolution 
as a merged object consists of only part of the RGB, but 
the full HB and AGB. 

In Fig. [121 we show a present-day, binned, theoretical 
HR diagram for our normal single-star population. Com- 
parison of Fig.[12]with the corresponding plot for merged 
objects (Fig. [7^,) indicates little difference between the 
two plots. This apparent similarity results from the use 
of our grid of normal single-star models to calculate the 
structure and evolution of the merged objects. This ap- 
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Total mass of single stars at present epoch (M Q ) 

Figure 11. Distribution of the total masses in a comparison 
model of normal present-day single stars constructed from the 
same grid of stellar evolution models used for the merged popu- 
lation. As in the corresponding plot for the merged population 
(Fig.Ef, the distribution is separated into three sub-populations: 
RGB stars (dashed curve), HB stars (solid curve) and AGB stars 
(dotted curve). 
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Surface temperature of single stars at present epoch (K) 

Figure 12. Two-dimensional distribution of the luminosities and 
effective temperatures in a comparison model of normal present- 
day single stars constructed from the same grid of stellar evolution 
models used for the merged population. The grey scale represents 
the total number of merged objects per two-dimensional bin. Tem- 
perature increases to the left on the horizontal axis to facilitate 
comparison with observational HR diagrams. 



proximation was necessary because of the lack of suitable 
stellar structure and evolution models for objects that 
merge during the CE phase. The degree of error intro- 
duced in our predicted luminosities and effective temper- 
atures by using non-rotating stellar models is not clear. 

4. DISCUSSION 

We have shown that a distinguishing characteristic fea- 
ture of a merged stellar population is the high rotational 
velocity in comparison with a single-star population. 
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Specifically, rapid rotation in the merged population is 
expected for a fraction of stars in their giant phase or 
horizontal-branch phase of evolution. Among observed 
RGB stars, rapid rotatorqj account for 1% to 6% of the 
red giant populatio n, depending on the specific obser- 
vational sample (see iCarlberg. Majewski. fc Arrasll2009l . 
and references therein). Fig. [2] shows that ~70% of the 
projected rotational velocities in our standard merged 
object population are greater than lOkms -1 . Combin- 
ing this with the number of RGB stars in our merged 
object and normal single star model populations from 
Table [21 and assuming a ZAMS binary fraction of 0.5, 
we estimate the fraction of rapid rotators in the RGB 
population due to CE mergers is ~ 3.4%. We note that 
this estimate neglects RGB stars in binaries whose orbits 
are wide enough to avoid mass transfer, since such bina- 
ries are excluded from our population of merged objects. 
Inclusion of these binaries increases the overall number 
of RGB stars, which would lower our estimate. How- 
ever, it is also true that tidal interactions between the 
RGB star and its lower-mass companion can still spin up 
the giant star in some of these detached binaries. Even 
though such a system retains its binarity, the presence 
of the companion may be masked by the much brighter 
RGB star and the system could be misidentified observa- 
tionally as a solitary RGB rapid rotator. At the moment, 
it is uncertain to what extent such misidentifications are 
present in the observational estimates above. We also 
point out that ICarlberg. Majewski. fc Ar ras (|2009f ) esti- 
mate that rapid rotators produced by the ingestion of 
planets by RGB stars, which we do not consider, can 
account for ~ 0.6% of the RGB population. 

The subsequent evolution of rapidly rotating RGB 
stars can lead to the production of rapidly rotating HB 
stars as a result of the contraction of stars from the tip 
of the RGB to the HB when helium is ignited in the 
core (see §[3]). Although we find theoretical evidence for 
the rapid rotation of stars in a clump on the red portion 
of the HB in our HR diagram (see Fig. [7]), observational 
evidence for rapid rotation among red HB field stars is 
lacking (e.g.,|Behr 2003; dc Mcdciros 200-1). However, in 
a sample of 45 HB field stars. iBehrl ()2003t ) found evidence 
for rapid rotation (v sin i > 30 km s™ 1 ) in several blue HB 
field stars (7500 K < T < 11,500 K). He suggests that the 
underlying distribution of actual rotational velocities in 
blue HB field stars may be bimodal, consisting of a slow 
population (wsini ~ 10 km s~ x ) and a fast population 
(i>sini ~ 30 km s _1 ), similar to the bimodal distribution 
found in many globular clusters. If a larger sample of HB 
field stars confirms the possibility of a bimodal distribu- 
tion in rotational velocities, the fast-rotating population 
may result from the mergers studied in this paper. As al- 
ready pointed out bv iPolitano et al.l (|2008f ). the merged 
population of very rapidly rotating HB stars may evolve 
under the action of enhanced stellar wind mass loss due 
to the reduction of surface gravity associated with cen- 
trifugal effects. As a result, some very rapidly rotating 
HB stars may lose a significant fraction of their envelope, 
which would cause both a reduction in the rotation rate 
of the star because of the associated angular momentum 
loss, and a blueward movement in the HR diagram. Such 

6 We follow the literature in defining rapid rotators as those 
RGB stars with vsvai > lOkms -1 . 



stars may contribute to t he rapid rot ation of blue HB 
stars in the field found by IBehrl (120031 ) and, if sufficient 
mass is lost, provide an ev olutionary channel for the for- 
mation of single sdB stars. IPolitano et al.1 (J2008I ) predict 
a spectrum of masses ranging from 0.32 to ~ 0.7 Mq with 
a strong peak between 0.47 and 0.54 Mq. This is consis- 
tent with the small number of observationally determined 
masses for single sdB stars, the majority of which have 
masses between 0.46 and 0.54 Mm, but can be as small 
as 0.39 Mq (PC 0911+456: lRandall et alJl2007ft . We fur- 
ther note that, although our models have solar mctallic- 
ity, if similar rapid rotation and subsequent mass loss on 
the HB is found for merger models with lower metallici- 
ties, this may be relevant to the exten ded HB morpholo- 
gies found in globular clusters (e.g., iMoni Biden et all 
l20t)8ft . 

Another class of stars that possibly may result from a 
merger during CE evolution are FK Comae stars. These 
stars are rapidly rotating gia nts with vsini ~ 1 00 km s ~ 1 
and spectral types G and K (jBopp fc Stencelll"l98lD . An 
origin involving the evolution of a rapidly rotating sin- 
gle MS star is not a viable one, since the progenitors 
of FK Comae stars would have rapid rotation and be 
chromospherically active, whereas none have been dis- 
covered. On the oth er hand, the co alescence of evolved 
contact binaries (see lWebbinklll976T) remains as a viable 
evolutionary scenario, and we suggest that the merger of 
non-corotating CE binaries may also contribute to this 
population. We note, however, that our current model 
population does not contain merged objects of the cor- 
rect spectral type, although this may be due to our use 
of non-rotating stellar models for these objects. Incor- 
poration of stellar models that include rotation would be 
especially important in determining the viability of the 
merger scenario for FK Comae stars. 

As a consequence of the merger process, it is likely 
that there exists a nonsphcrical distribution of circum- 
stellar matter surrounding the remnant, either resulting 
from matter which is not accreted nor ejected from the 
system during merger or from an enhanced mass-loss 
phase from the equatorial region of the rapidly rotat- 
ing star. In the former case, there may exist a short 
phase (depending on the mass in the circumstellar enve- 
lope) where the interaction of the stellar wind with this 
disk structure can lead to the development of a bipolar 
outflow (where outflow in the equatorial region is im- 
peded by the presence of the circumstellar disk). Fur- 
thermore, it is likely that the mass-loss history of the 
remnant following the merger event would be variable, 
and that the circumstellar matter distribution would re- 
flect its integrated history since the merger event. For 
example, such a temporal history may be responsible for 
the formation of a torus-like structure embedded within a 
continuous outflow, the existence of which has been pos- 
tulated for the shaping of planetary nebula e in the gen- 
eralized interacting- stellar-winds model (see lKwoklll982l : 
iKahn fc~W est 1985]), and which has been inferred from 
HST observatio ns of the proto-p lanetary nebula IRAS 
17106-3046 (see lKwoket alll200(T . 

Such structures are of interest, especially for AGB 
stars, since observational evidence of axisymmetric 
circumstellar envelopes in po st-AGB stars and pre - 
planctary nebulae exists (e.g., iSahai fc Traugerlll998l) . 
The mass loss in AGB stars is generally spherically sym- 
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metric (e.g., iBowers et alj|1983t ISahai fc Bieglnel Il993t 
Olofsson 2004) , although observational evidence of asym- 
metrical structure exists in the AGB star V Hya (see 
iHirano et al.l 12004 ). Given this state of affairs, the evo- 
lutionary state of stars leading to the transition from a 
spherical circumstellar envelope to a nonspherical mor- 
phology is of great interest both theoretically and obser- 
vationally. The results of our study suggest that the 
frequency of occurrence of nonspherical circumstellar- 
envelope structures resulting from the merger of the com- 
ponents of a binary system during CE evolution is not 
high (at best ~ 10%, assuming all merged AGB stars 
contain asymmetric envelopes, although it is not clear 
whether mergers between AGB stars and planets would 
lead to asymmetric envelope structures as well). Hence, 
if the nonspherical structures are not uncommon and a 
result of binary star interactions, then they may reflect 
the population of det ached systems in wide binaries {e.g., 
iHuggins et al.l l2009f ) or those remnant binary systems 
that have survived the CE phase resulting in short or- 
bital periods. An extension of our population synthesis 
study to examine the orbital properties of the detached 
systems and the effect of gravitational focusing on the 
shape of the AGB envelopes will be the subject of a fu- 
ture investigation. 

The effect of rotation in the deeper interior layers of the 
merger remnant may affect the composition of matter in 
the surface layers if rotationally induced mixing of mat- 
ter develops. In such cases, the elemental abundances 
in the surface layers may show anomalies compared to 
those of normal single stars, as matter may be mixed 
from the outer radiative regions of the hydrogen burning 
shell to the base of the convection zone and, hence, to 
the surface. The changes in the abundances, however, 
sensitively depend on details regarding the magnitude of 
the diffusion coefficient and the d epth of the mixing (see 
for example. iPalacios et al J 120061 ). 

Of particular interest with regard to the nucleosyn- 
thetic anomalies at stellar surfaces are the Li-rich stars. 
These stars are located near the red-giant bump region at 
luminosities of ~ 60 L©, and i n the early AGB phase at 
luminosities of ~ 600 L © (see IReddv fe Lambert! 120051 : 
Kumar fc Reddvl 12009ft . The stars near the red-giant 
bump may have masses of the order of ~ 1 — 2.4 M Q , 
based on comparison of their luminosities and effective 
temperatures with red-giant evolutionary tracks. A com- 
prehensive discussion of the theoretical models (includ- 
ing the merger of a red giant with a Jupiter mass type 
planet), involving the conversion of 3 He to 7 Be, the for- 
mation of 7 Li via electron capture (jCameron fc Fowlerl 
|1971| ), and mixing associated with an enhanced diffu- 
sion coeffici ent (presumably a resul t of ra pid rotation), 
is given in [Pcnisscnk ov fc Herwid (|2004l) . If the ef- 
fectiveness of the rotation-induced mixing is related to 
the location of the bottom of the surface convection 
zone relative to the hydrogen burning shell, then such 
mixing may take place in stars above the base of the 
giant branch, as the mass lying between the base of 
the convection zone and the outer edge of the hydro- 
gen burning shell decreases with increasing evolution- 
ary state. Given that the proportion of Li-rich giants 
among rapid rotators (vsmi >8kms _1 ) is ~ 50% (see 
Dcnissenkov & Pinsonneault 2008), which is consider- 



ably greater than the ~ 2% of Li-rich stars among slowly 
rotating stars (wsini < lkms -1 ), it is suggestive that 
rapid rotation may play an important role in forming a 
significant fraction of the Li-rich giant population. 

An additional observation suggesting a possible con- 
nection between Li-rich stars and the merged remnants 
envisioned here is the relation between the Li abundance 
and the far-infrared e xcess in these stars. In particu- 
lar, Ide la Reza et al.l (|1996l ) suggest that an enhanced 
mass loss event may accompany the Li-enrichment phase. 
The mass loss resulting from the phase is expected to be 
significantly greater (about 2 orders of magnitude) than 
that for normal Li-poor giants. It should be noted that 
the far-infrared emission may have its origin in a disk-like 
configuration resulting from the merging process rather 
than an expanding dust shell. 

5. SUMMARY AND CONCLUSIONS 

The population of binary star systems consisting of an 
RGB or AGB star with a lower-mass companion that 
merge during CE evolution has been explored via a pop- 
ulation synthesis technique. Using a suite of 116 evolu- 
tion tracks for single stars and including the rotational 
evolution of the primary components, we have evolved 
10 7 ZAMS binary systems. The characteristic properties 
of the merger remnants have been examined as a func- 
tion of uncertain population synthesis input parameters, 
such as the initial mass-ratio distribution of the binary 
population and the efficiency for mass ejection in the CE 
phase, as well as the critical rotational velocity at which 
enhanced mass loss occurs for the merged object. Using 
simplified prescriptions for the binary interaction during 
the CE stage and for the rotational and nuclear evolu- 
tion of the merger remnant, we find that the present- 
day merged population is primarily distinguished by its 
high rotation speed. For the total merged population, 
the median projected rotational velocity is 16.2 km s -1 , 
in contrast to 2.3kms~ 1 for a theoretical population of 
normal single stars calculated using the same stellar mod- 
els and initial mass function. This relative difference is 
also reflected in individual sub-populations of merged ob- 
jects (see below), whose peak rotational velocities are an 
order of magnitude higher than those in corresponding 
sub-populations of normal single stars. 

In our present-day model population of merged ob- 
jects, we include only those objects that are on RGB, HB 
or AGB. We find that these objects constitute between 
0.24% and 0.33% of the initial population of 10 7 ZAMS 
binaries, depending upon the values chosen for the input 
parameters. In addition, we find that HB stars are the 
most populous (57%), followed by RGB stars (37%) and 
then AGB stars (6%). The luminosities of these sub- 
groups overlap and lie in the range of ~ 10 — 100 L@, 
with the majority of the present-day merged population 
consisting of HB stars with luminosities of ~ 60 L Q , vir- 
tually indistinguishable from their present-day, normal 
single-star counterparts with the exception of their high 
rotation speed. The masses characteristic of the total 
population in our standard model are < 2 M Q , with me- 
dian masses for RGB (1.20 M Q ), HB (1.35 M ) and AGB 
(1.34 Mq) stars that do not significantly differ from those 
for normal single stars on the RGB, but are clearly lower 
than those for single stars on the HB and AGB. The me- 
dian mass of the total population is about 0.05 Mq higher 
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for our standard merger population than for a normal 
single-star population, which may seem to contradict the 
results per sub-population, but can be explained by the 
different RGB versus HB ratios between merged objects 
and normal single stars (see § I3.3|) . These quantitative 
features of the present-day merged population arc found 
to be relatively insensitive to the population synthesis 
input parameters. 

The results of our population synthesis study have been 
discussed in terms of possible observational counterparts 
either directly involving the high rotational velocity of 
the component or indirectly, via the effect of rotation 
on the amount and distribution of circumstellar matter. 
Additional probes potentially indicating the occurrence 
of such interactions could be provided by compositional 
signatures on the stellar surface. As a specific application 
of our CE merger model, we have considered the possi- 
bility that rotationally enhanced mass loss may be suf- 
ficient to significantly deplete low- or intermediate-mass 
hydrogen-rich envelopes of rapidly rotating HB stars to 
the extent that their effective temperatures are increased 
sufficiently to produce sdB stars. Such a merger sce- 
nario can provide an additional evol utionary channel 
for the formation of single sdB stars (jSoker fc Harpazl 
120001 120071: iPolitano et all 12008ft . complementary to the 
double hel i um- W D merger channel proposed earlier by 
lHanet al.1 (12003ft . 

Indirect effects associated with the merger and rapid 
rotation, coupled with the evolution of the merged object 
to its later phases, can lead to deviations from spheri- 
cal symmetry in the distribution of circumstellar mat- 
ter. Such structures can be potentially prob ed via their 
emission in the radio (e.q.. lHirano et aLll2004ft or infrared 
wavel ength region, or via the s cattering of optical light 
(e.g., iMauron fc Huggiiisl 12006ft . These structures may 
result from inherent asymmetries present in the merger 
process itself, or result after the merger has already oc- 
curred, reflecting asymmetric mass loss from the stellar 
surface or the interaction of a spherically symmetric stel- 
lar wind with a nonsphcrical distribution of circumstellar 
matter. Additional observational effects may be a conse- 
quence of rotation in the deep interior, possibly resulting 
from rotationally induced compositional mixing leading 
to the production of nucleosynthetic anomalies at the 
stellar surface. In this regard, the merger scenario may 
provide an additional evolutionary channel for producing 
Li-rich RGB stars. 

To gain insight into the outcome of the merger sce- 
nario, our investigation has been primarily exploratory 
in nature. As a result, the study has been based on 
a number of simplifying assumptions. Foremost among 
them is the treatment of the structure and evolution of a 
merged object as equivalent to a single, non-rotating star 
of the same total mass and core mass. The binary inter- 
action involving the evolution into and during the CE 
phase is also uncertain, as the efficiency of mass ejection 
has been assumed to be a constant, independent of evolu- 
tionary phase. Mass loss prior to and after the CE phase 
has been based on rates inferred for normal single stars, 
which are likely to have much smaller rotational veloci- 
ties. Finally, it is likely that the mass loss via a stellar 
wind, and hence angular momentum loss from the sys- 
tem, would be affected by the spin-orbit tidal coupling 
and magnetic braking. 



To improve on the numerical modeling, future theoret- 
ical work is indicated, especially in including neglected 
effects. For example, the use of a grid of stellar mod- 
els resulting from the merger process (see for example, 
iGlebbeek fc Pola 12008ft would be particularly illuminat- 
ing, not only in describing the structure and properties of 
the merged remnant, but also in determining their evolu- 
tionary characteristics as a function of age. This would 
be most useful for our large population of merged ob- 
jects on the RGB and HB. Our work also provides valu- 
able guidance by better defining the region of parameter 
space where detailed multi-dimensional hydrodynamical 
calculations are highly desirable. We find that mergers 
during CE evolution likely happen in binaries character- 
ized by primaries with masses of ~ 1 — 2 Mq that fill their 
Roche lobes on the lower half of the RGB at orbital pe- 
riods between ~ 1 and 30 days. Our study also points to 
the phases of evolution where hydrodynamical or magnc- 
tohydrodynamical calculations are needed to explore the 
influence of rotation and/or magnetic fields in the mix- 
ing of matter and angular momentum from the the deep 
interior (e.g., regions that have undergone nuclear pro- 
cessing) to the stellar surface. Incorporation of mixing 
prescriptions, based on the above studies, into evolution- 
ary calculations of a merged remnant could provide quan- 
titative estimates of the influence of rotation on the sur- 
face composition. Of special interest would be numerical 
modeling of the rotating remnant during a helium flash, 
as it may provide a compositional probe of the possible 
mixing of matter between the helium core and hydrogen- 
rich envelope for stars that are now on the HB. Finally, 
an extension of our models to examine the properties of 
binaries that merge in the MS phase of evolution is of in- 
terest, as they may have applications to highly luminous 
transient events, perhaps similar to that observed in V 
838 Mon (see for example. ISoker fe Tvlendal 12006ft . 

In addition to the fundamental theoretical work, fu- 
ture observational programs aimed at detecting the ef- 
fects of rapid rotation using ground-based optical inter- 
ferometers could be fruitful in the determination of the 
shapes of giant stars, where tem perature variation s of 
the surface may be probed (e.g., iZhao et al.ll2009ft . A 
program of highly accurate photometric studies can be 
considered, enabling one to detect differences in lumi- 
nosity and colors of rapidly rotating stars in comparison 
to their single-star counterparts. For example, observa- 
tional evidence exists that magnetically-active, rapidly 
rotating, low-mass MS stars in eclipsing binaries are 
characterized by effective temperatures that are lower 
than, radii that arc larger than, and luminosities that are 
approximately equal to those predicted from correspond- 
ing theoretical ste llar models of the same mass (e.g., 
iTorres et al.l 12009ft . The same stellar models satisfac- 
torily match observations of similar-mass s tars that are 
rotati n g more slowly in wide binary orbits. | Torres et al.l 
(|20"06ft . lL6pez-Moralesl ([2007ft . iMorales et all (|2008ft and 
others have shown that the source of these differences 
is the magnetic activity in the rapidly rotating stars. 
Whether the radii, effective temperatures and luminosi- 
ties in rapidly rotating giant or HB stars would be sim- 
ilarly affected is an interesting and unresolved question. 
Finally, studies on the asymmetries in the circumstel- 
lar envelope of stars in the AGB phase with upcoming 
facilities (e.g., ALMA) will be especially useful in provid- 
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ing constraints on the stellar evolutionary phases where 
asymmetries develop. Such studies will be very impor- 
tant for distinguishing the contributions of the merged 
population from the existing binary population, thereby 
potentially providing further constraints on uncertain 
population synthesis input parameters and possible shap- 
ing mechanisms for the origin of the asymmetries seen to 
be prevalent in the post-AGB and proto-planetary neb- 
ula phase. 
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